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Abstract 

The nonparametric test for change-point detection proposed by Gombay and 
Horvath is revisited and extended in the broader setting of empirical process theory. 
The resulting testing procedure for potentially multivariate observations is based on 
a sequential generalization of the functional multiplier central limit theorem and on 
modifications of Gombay and Horvath's seminal approach that appears to improve 
the finite-sample behavior of the tests. A large number of candidate test statistics 
based on processes indexed by lower-left orthants and half-spaces are considered and 
their performance is studied through extensive Monte Carlo experiments involving 
univariate, bivariate and trivariate data sets. Finally, practical recommendations 
are provided and the tests are illustrated on trivariate hydrological data. 

Keywords: half-spaces; lower-left orthants; multiplier central limit theorem; multi- 
variate independent observations; partial-sum process. 
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1 Introduction 



Let Xi, . . . ,Xn be a sequence of independent (i-dimensional random vectors for some 
fixed integer d ^ 1. The aim of this work is to study, both theoretically and empirically, 
nonparametric tests for the detection of a change-point in the sequence Xi, . . . , X^- The 
corresponding null hypothesis is 

Hq : 3 Pq such that Xi, . . . , X„ have law Pq. (1) 

As frequently done, the behavior of the derived tests will be investigated under the 
alternative hypothesis of a single change-point: 

Hi : 3 distinct Pi and P2, and k* s [1, . . . ,n — 1} such that 

Xi, . . . , Xk* have law Pi and . . . , X^ have law F2- (2) 



There exists an abundant literature on nonparametric tests for change-point detec- 
tion. We shall not review here procedures designed for serially dependent observations. 
The approaches proposed for sequences of independent observations differ, on one hand, 
according to the test statistic, and on the other hand, according to the resampling tech- 
nique used to compute an approximate p-value for the test statistic. In terms of the 
test statistic, two frequently encountered classes of approaches are those based on U- 
statistics (see e.g. ICsorgo and Horvathl . llQSSt iFergerl . Il994t iGombav and Horvathl. l2002l: 



1999 



Horyath and Huskoval. 120051) a nd those based on empirical c.d.f.s (see e.g. lGombav and Horvathl . 



Horvath and Shaol . l2007h . As far as the resampling technique is concerned, one finds 



approaches based on permut a tions of the original sequence (s ee e.g. lAntoch and Huskova 



200 ll : iHorvath and Huskoval . l2005l : iHorvath and Shad . l2007l ) and ap proaches that use a 



Gombav and Horvath . 



weighted b ootstrap based on multiplier central limit theorems (see e.g 
I999I . l2002l ) . For a broader pr esentation of the field of change-p oint analysis, we refer the 
reader to the monographs by iBrodskv and Darkhovskvl f 19931 ) and ICsorgo and Horvath 



(119971) 



In this paper, we revisit and extend the approach proposed by iGombay and Horvath 



( I1999I ) based on the test statistic 



T 



max 

l<fc<n-l 



k{n — k) 



n 



3/2 



sup Ffc(x) 



X] 



where 



-| fc 1 " 

Fu{x) = - 2 < ^) and FU{x) = -— ^ J] 1(X, < 



i=l 



X) 



X e 



j=fe+i 



are the empirical c.d.f.s compu ted from Xi, . . . ,Xk and X^+i, . . . , respectively (see 



also lCsorgo and Horvathl . Il997l . Section 2.6). From a theoretical perspective , we work in 



(2000 


) and Kosorok 


( 


2008) 



of functions (in the sense of empirical proce ss theory - see Section [2?T|) . we first extend 



the multiplier central limit theorem (see e.g. iKosorokl . l2008l Theorem 10.1 and Corollary 
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10.3) to the sequential settin g. This allows us to o btain interesting generalizations of 
Theorems 2.1, 2.2 and 2.3 of iGombav and HorvathI (119991 ). In particular, we propose a 
slightly different multiplier process that appears to lead to better behaved tests in the 
case of moderate sample size. From a more practical perspective, we consider a large 
number of candidate test statistics based on processes indexed by lower-left orthants and 
by half-spaces, and we study the finite-sample performance of the corresponding tests 
through extensive Monte Carlo experiments involving univariate, bivariate and trivariate 
data sets. As we shall see, in the multivariate case, the tests based on processes indexed 
by half-spaces appear to be substantially more powerful than more classical tests based 
on multivariate empirical c.d.f.s (i.e., based on processes indexed by lower-left orthants). 

The paper is organized as follows. In the second section, we state the theoretical 
results at the root of the studied class of tests in the broad setting of empirical process 
theory. The third section is devoted to an application of the theorems of Section [2] to the 
derivation of nonparametric tests for change-point detection for two classes of functions 
which are the collection of indicator functions of lower-left orthants and the collection 
of indicator functions of half-spaces. The results of large-scale Monte-Carlo experiments 
comparing the finite-sample behavior of the tests are partially reported in the fourth 
section. The last section contains practical recommendations and presents an application 
of the studied tests to trivariate hydrological data. All the proofs are relegated to the 
appendices. 

Note finally that the code of all the tests studied in this work will be documented and 
released as an R package whose tentative name is npcp. 



2 Theoretical results for change-point detection 
2.1 Notation and setting 

All the random variables used in this work are defined with respect to the underlying 
probability space {fl, Q, P) and the outer probability measure corresponding to P is de- 
noted by P*. 

Let Xi, . . . , Xn be i.i.d. d- dimensional random vectors with law P, and let J-" be 
a class of measurable functions from M'^ to M. The empirical measure is defined to 
be P„ = ^^^i6xi, where Sx is the measure that assigns a mass of 1 at x and zero 
elsewhere. For f ^ J^, Pn/ denotes the expectation of / under Fn, and Pf the expectation 
under P, i.e., 



1 " r 

Pn/ = -y]/TO and Pf=\fdP 



The empirical process evaluated at / is then defined as G„,/ = ■y/n(Pn.f — Pf)- 

Saying that J-" is P-Donsker means that the sequence of processes : / e 

converges weakly to a P-Brownian bridge {Gp/ : / e J^} in the space i'^{J^) of bounded 
fu nctions from to R equipped wit h the uniform metric in the sense of Definition 1.3.3 



of Ivan der Vaart and Wellnerl (120001 ). Following usual notational conventions, this weak 



convergence will simply be denoted by G„ -w Gp in £°°(J-'). Furthermore, we say that 
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Fe : M is an envelope for T if -Fg is measurable and < Ff.[x) for every j ^ T 

and X G W^. 

The advantage of working in this general framework is that the forthcoming results 
remain valid for many P-Donsker classes T . By taking T to be the class of indicator 
functions of lower-left orthants in M'^, i.e., T = {y ^ \{y x) : x £ R } with M = 



Mu {—00, oo}, one recovers the setting studied in lCsorgo and HorvathI ( 1l997l . Section 2.6) 
and based on empirical cumulative distribution functions (c.d.f.s). Although this is a 
natural choice for J-", many other choices might be of interest in pra ctice such as the class 
of indicator functions of closed balls, rectangles or half-spaces (see iRomand . Il988l . for a 
related discussion regarding the choice of J-"). 



2.2 A multiplier central limit theorem for the sequential empir- 
ical process 

The sequential empirical process is defined as 

Z„(s, /) = ^ y,{f\X,) - Pf) = VK{^)GynsiL s e [0, 1], / 6 ^, 

where A„(s) = [nsj/n and with the convention that Pq/ = for all j ^ T . 

According to Theorem 2.12.1 of Ivan der Vaart and Wellnerl ( 120001 ) . T being P-Donsker 
is equivalent to Z„ Zp in ^^([0, 1] x J^), where Zp is a. tight centered mean-zero Gaus- 
sian process with covariance function 

coy{Zp{s,f),Zp{t,g)} = {s a t){Pfg-PfPg) 

known P-Kiefer-Miiller process. 

Given i.i.d. random variables ^i, . . . , ^„ with mean and variance 1, satisfying 5[j°{P(|^i| > 
a;)}^/^dx < 00, and independent of the random sample Xi, . . . , X„, we define the following 
multiplier yeision of Zn- 

Z,(5, f) = ^y ^dfiX^) -Pf}, s e [0, 1], / E 

Notice that the empirical process Z„ depends on the unknown map / t-^- Pf and therefore 
cannot be computed. With applications in mind, we define two versions of Z„ (depending 
on how / 1-^ Pf is estimated) that can be fully computed. For any s e [0, 1], / g J^, let 

[ns\ [ns\ 



4=1 ^ J=l 



where ^[^^j = [ns\ ^ YjI^i = by convention, and let 

[ns\ 

z„(s,/) = ^Ve.{/ra-Pn/}. 



n . ^ 
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The following result is then a partial extension of the multiplier central limit theorem 
(see e.g. iKosorokl . l2008l . Theorem 10.1 and Corollary 10.3) to the sequential setting. 



Theorem 1. Let T he a P-Donsker class with measurable envelope such that PF^ < 
00. Then, (Z„, Z.„, Z„, Z„) -w (Zp, Zp, Zp, Zp) in {£°°([0,1] x J")}^, where Zp is an 
independent copy ofZp. 

Theorem [T] suggests the following interpretation: when n is large, Z„ can be regarded 
as "almost" an independent copy of Z„, while Z„ and Z„ can be regarded as computable 
copies of Z„. As we shall see, this interpretation is at the root of the resampling technique 
considered in Section [31 

Although each of Z„ and Z„ could be regarded as "almost" an independent copy 
of Zjj, their behavior for moderate n might differ quite substantially. In Section HJ we 
empirically investigate which of Z„ or Z„ leads to tests for change-point detection with 
the best finite-sample properties. 



2.3 Application to change-point detection 

Recall that the null and alternative hypotheses under consideration are given in ([1]) 
and (EJ, respectively. 

Let J-" be a class of measurable functions. In o rder to test the aforementio ned hypothe- 
ses, we extend the approach studied in detail by lCsorgo and HorvathI (Il997l . Section 2.6) 
and compare, for all A; e {1, . . . , n — 1}, 

^kf = lj]f{X.) and = -L^ f(X,), f e F. 



i=l 



i=k+l 



Analogous to ICsorgo and HorvathI (Il997l . Section 2.6), we define the process 

/) = ^\n{s) {1 - A„(S)} (P^„,j/ - P:_^„,j/) , 5 G [0, 1] , / G J", 

where A„(s) = [nsj/n and with the convention that Pq/ = and Pq/ = for all f e F. 
Notice that, under the null hypothesis, for any s s [0, 1] and f e F, we have 

/) = {1 - Xn{s)} Zn{s, f) A„(s){Z„(l, /) - Z„(s, /)} = Z„(s, /) - A„(s)Z„(l, /). 

(3) 

With resampling in mind, we define two multiplier versions of D„ based on the mul- 
tiplier versions of Z„ defined in the previous subsection. For any s e [0, 1] and f e F, 
let 

/) = {1 - A„(s)}Z„(s, /) - A„(s){Z„(l, /) - Z„(s, /)} = Z„(s, /) - A„(s)Z„(l, /), 
and, following iGombay and Horvath (Il999f ). let 

D„(s, /) = {1 - A„(s)}Z„(s, /) - Xn{s)Kis, /), (4) 
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where 



KisJ) = ^ y. {^i-C-[ns\)fi^i) with C-[ns\ = j T J] Ci, (5) 

^Jn f—^^ ^ ^ ^ ^ n — \ ns\ f—^^ 

* i=[nsj + l ^ ^ i=\ns\-\-\ 

and ^0 = by convention. Notice that the process Z* defined above is, up to a small 
error term vanishing as tt, ^ oo, the version of the process (s, /) ^ Zn(l ^ ■§) /) computed 
from the "reversed" sequence X„_i), . . . , (6, Xi). 



The following two results extend Theorems 2.1, 2.2 and 2.3 of iGombay and Horvath 



( I1999I ) and suggest, for large n and under the null hypothesis, to interpret each of D„ and 



D„ as an "almost" independent copy of D„. 

Theorem 2. Assume that Hq holds and that T is P^-Donsker with measurable envelope 
Fe such that PqF^ < 00. Then, (D„, D„, D„) ^ (Dp^, D^^, D'pJ m {£^([0, 1] x T)}^, where 
Dpq(s, /) = Zpq(s, /) — sZpq(1, /), s s [0,1], j ^ T , and Dp^ is an independent copy of 

As we continue, for any J : J-" M, supjgjr \ Jf\ will be denoted by || Also, for any 
sequence of maps Yi,Y2, . . . , each from Q to M, we say that the sequence Yn is bounded in 
outer probability if, for any e > 0, there exists M > such that sup^gj^ P* > M) < e. 

Theorem 3. Assume that Hi holds with k* = [nt\ for some t e (0, 1) and that T is P\ 
and P2-Donsker with measurable envelope such that PiF^ < 00 and P2F^ < 00. Then, 

(t) sup,g[o,i] \\n-"^^n{s, f) - /)||^ ^ 0, 

where Kt{s, f) = {PJ - P2f){s a t){l - (s v t)], 

(a) supsg|-Q i] ||D„(s, is bounded in outer probability, 
(Hi) i^n converges weakly in £^([0, 1] x J^). 

The previous result will be used in the next section to show that various tests for 
change-point detection based on D„ will tend to reject Hq under Hi as n increases. 

3 Tests for change-point detection a la Gombay and 
Horvath 

The aim of this section is to use the results of the previ ous section to derive t e sts fo r 



change-point detection in the spirit of those proposed by iGombay and HorvathI (jl999f ) . 



Among the many possible choices for J-", we consider the following two: 
(CI) the collection O of indicator functions of lower-left orthants in M"', where 

O = {/,(?/) = 1(2/ <x):x el'}; 
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(C2) the collection T-L of indicator functions of half-spaces in W^, where 

^ = [faM = Ua'^y ^ 6) : a E 5,, 6 e 1} 
and Sd is the subset of M'^ composed of vectors with Euclidean norm one. 

It is well-known that lower-left orthants and half-spaces are Vapnik- Chervonenkis classes 

of sets. Consequently , O and T-L are P-Donsker for any law P (see e.g. Ivan der Vaart and Wellnerl . 



2nnnl : lKosorokl . l2nn8h . 



As we continue, in the case of choice (CI), for any s e [0, 1] and fx e O, ^n{s, fx), 
D„(s, fx) and ]D)„(s, fx) will simply be denoted by D„(s, x), D„(s, x) and D„(s, x), respec- 
tively. Similarly, in the case of choice (C2), for any s e [0,1] and fa^b ^ T^, ^n{s,fa,b), 
D„(s, fafi) and B>n{s, fa,h) will be denoted by D„(s, a, 6), D„(s, a, h) and ©^(s, a, 6), respec- 
tively. 

In the framework under consideration, a change in the sequence Xi, ■ ■ ■ ■, X^ can occur 
at any point e {1, . . . , n — 1}. A test for change-point detection could therefore be 
obtained by first defining a test statistic for any possible change-point /cE{l,...,n — 1}, 
and then by combining the resulting n — 1 statistics into a global statistic using some 
function from ip : M"^-*^ M. 

In the case of choice (CI), two natural possibilities for the n — 1 change-point statistics 
are respectively 



where Fn{x) = fnfx, a; e M , is the empirical c.d.f. computed from Xi, . . . , X„, and 



Tn^k = sup 











= max 











k e {1, . . . ,n — 1}. 



Two natural choices for the function -0 are the maximum and the arithmetic mean which 
leads to the following four global statistics: 



Sn,v = max Sn,k = sup 
r„.v = max Tr 



n,v — j.j.±ci^ -tn,fc — sup sup 

n-1 



'„ {s,x)}'^ dFn{x), 
I (s, x) I , 



1 n-l r.1 r. 

S^^+ = -ySn,k= \ {D4s,x)}'dF„(x)ds, 
k^i Jo 

^ n-1 ^1 
^n,+ = ~ / I Tn k = 



sup 



K {s,x)\ds. 



Note that r„ v is the global statistic used in iGombay and HorvathI fjl999l ). 
In the case of choice (C2), for any k e {l,...,n — 1}, we first define 











= sup max 






ae5+ 





where = {a e Sd '■ ai > 0} and, for any a e , Fa,n is the (univariate) empirical c.d.f. 
computed from the projected sample a^Xi, . . . , a^Xn, and 

K,fc = sup sup 

aeS+ beR 

As in the case of choice (CI), this leads to four global statistics depending on whether the 
change-point statistics are combined using the maximum or the arithmetic mean, i.e., 

Un,v = max Un,k = sup {Dn{s,a,b)f dFa,nib)da, 

l^fcs;n-l ^g[o_;^] Jg+ Jjj 

K,v = max K,fc = sup sup sup (s,a, 6)1 , 

Un,+ = -y,Un,k=\\ f {Dn (s, tt, 6)}' dF,,„(6)dads, 
" fc=i ^0 Js+ Jk 

K,+ = - Y! K,fc = sup sup |D„ (s, a, 6)1 ds. 

In our Monte Carlo experiments, the integral and the supremum over a e in the 
definitions of Un,k and Vn,k, respectively, were approximated numerically based on a 
uniform discretization of into m points. The implementation of the tests based on Un,k 
and Vn^k is discussed in more detail in Appendix [Bl Notice finally that the change-point 
statistics Sn,k and Un^k (resp. T„ fc and Vn^k) coincide when d = 1 since = {1}. 

Let us now explain how approximate p-values for these statistics can be computed 
using the multiplier processes D„ and D„. For the sake of brevity, we present the approach 
and state the key results only for 5'„,v as the cases of the other seven global statistics are 
similar. 

Let be a large integer and let i e {1, . . . , n}, j g {1, . . . , N}, be i.i.d. random 
variables with mean and variance 1 satisfying 5(^{P(|e?| > x)y/Mx < 00, and indepen- 
dent of Xi, . . . , Xn. Also, for any j e {1, . . . , A}, let (resp. Dn '') denote the version 
of D„ (resp. D„) computed from ^[^\ . . . , Moreover, for any j e {1, . . . , A^}, let 

SL%= sup r {D(f)(s,x)}'dF„(x) and = g^p f {©(f) (s, a;)}' dF„(x). 

sg[0,1] JR"^ ^ ' se[0,l] JlRd ^ ^ 

The following result is then essentially a corollary of Theorem [2l 
Proposition 1. Under Hq, 

q{l) q{N) q{l) q{N)\ (q q{l) q{N) q{l) q{N)\ 

'-^71, V ) *-^ri, V ' ■ ■ ■ ' n, V ' n, v ' " ■ ■ ' '~^n, v I V^ v)*-'v''''' v )'-^v''''' v / 

in [0, oo)(^^+^\, where 



Sy = sup {Dp(,(s,x)} dFo{x) 

se[0,l] jRd- 

is the weak limit of Sn,v , Fq is the c.d.f. corresponding to Pq, and . . . , s[^^ are 
independent copies of . 
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The previous proposition suggests interpreting the Sn^v (resp. the Sn^v) under the null 
hypothesis as "almost" independent copies of Sn,v and thus computing an approximate 
p- value for S'^^v as 

]72]l(^SU^".v) or as ^2]l(^?v^^n,v). (6) 

Proposition 2. Assume that Hi holds with k* = [nt\ for some t e (0, 1). Then, Sn,v ^ 
+00 while, for any j s {1, . . . , N}, Sn^v and Sn^v are bounded in outer probability. 

A consequence of the previous proposition is that, under Hi, the approximate values 
for Sn,v will tend to zero in outer probability. As mentioned earlier, results analogous to 
Propositions [1] and |2] can be obtained for Sn,+, T^^y, t^n,v5 Un,+ , Ki,v and 



4 Monte Carlo experiments 



Large-scale Monte Carlo experiments were carried out in order to study the finite-sample 
performance of the tests defined in the previous section. Let Qn be one of S'„_v, Sn,+ , 
Tn^v, Tn^+, Un^^, Un,+, Ki,v and Vn,+. In the rest of the paper, the test based on Qn 
will be referred to as the test based on Qn (resp. Qn) when its approximate p- value is 
computed using the multiplier processes B>n^ (resp. ©n^). 

To compare the power of the aforementioned tests, several univariate and multivariate 
scenarios were considered and 1000 samples of size n were generated under each scenario 
for n E {50, 100,200}. Recall that the c.d.f. corresponding to Pq in Hq defined in ([1]) is 
denoted by Fq. Similarly, the distinct c.d.f.s corresponding to Pi and P2 in Hi defined 
in ([2]) will be denoted by Fi and F2, respectively, as we continue. In all scenarios, the 
multipliers appearing in the Dn^ and ©n^ were taken from the standard normal distri- 
bution. All approximate p-values were computed from A^ = 1000 multiplier realizations 
and the tests were carried out at the 5% level of significance. 



[Table 1 about here.] 



From the previous section, it is easy to verify that, in the univariate case, the change- 
point statistics Un,k (resp. Ki,fc) coincide with the Sn,k (resp. T„,fc) since Si = {1}, and 
that the Sn,k and the Tn,k are solely based on ranks. From the latter fact, it follows 
that, under Hq and the assumption of continuity of Fq, independent realizations of each 
of the four global statistics based on the Sn,k or the Tn,k can be obtained by computing 
these global statistics from independent samples of size n generated from the standard 
uniform distribution. A natural way of computing an approximate p- value for each of the 
four global statistics then consists of proceeding analogously to ([6]) using A^ independent 
realizations. As we continue, the resulting four univariate tests will be referred to as the 
tests based on S*^, 3*+, T* ^ and T^ + . 

Table [T] gives rejection percentages of Hq in dimension one for all the aforementioned 
versions of the tests when data are generated under Hq. As can be seen, the tests 
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whose approximate p-value is computed using the processes seem to be too hberal 
(at least for n s {50,100}), and more hberal than their version computed from the 
\ Nevertheless, as expected, the empirical levels of the multiplier tests improve as n 
increases in the sense that they become closer to the 5% nominal level. Note that the 
tests based on S"* y, 5* T^^ and T*^ provide a sort of benchmark as, by construction, 
they should hold their level well for any n provided is taken sufficiently large. 



[Table 2 about here.] 



Rejection percentages of Hq in the univariate case when data are generated under 
Hi are reported in Table El Three scenarios are considered: Fi and F2 are the c.d.f.s 
of the A^(0, 1) and the iV(0.5, 1) distributions, respectively; Fi and F2 are the c.d.f.s of 
the A^(0, 1) and the A^(0,2) distributions, respectively; and Fi and F2 are the c.d.f.s of 
the exponential E[l) and £'(0.5) distributions, respectively. Notice that the parameter t 
taking its values in {0.1,0.25,0.5} determines the position of the change-point in Hi as 
k* = [nt\. As one can see, the tests based on the processes ©n^ are consistently slightly 
more powerful than those based on the 'Dn\ while the rejections rates of the latter are, 
overall, comparable to those of the tests based on simulation from the standard uniform 
distribution. This merely appears to be due to the fact that the tests based on the 
processes Dn are slightly too liberal for the sample sizes under consideration. Notice 
that the differences in power decrease as n increases, as expected. Among the tests based 
on the &J\ the one based on Sn,v seems, overall, to be the best choice for detecting 
changes in mean, while the test based on T„^+ seems, overall, to be the best choice for 
detecting changes in variance. The former seems also to be the best choice, overall, when 
data are generated under the third scenario involving exponential distributions. If one is 
willing to make continuity assumptions on the underlying distributions, the tests based 
on S* y and T* _,_ are equivalently good candidates. Clearly, there exists more powerful 
test for change-poi nt detection if it is assurned that only a change in mean or variance 



can occur (see e.g. iBrodsky and Darkhovskyl . Il993l : ICsorgo and Horvathl . 119971 ). 



In dimension two and three, we considered multivariate distributions constructed from 



Sklarl fll959l )'s representation theorem. The latter result states that any multivariate c.d.f. 
F : [0, 1] whose marg inal c.d.f.s FW,...,FM are continuous can be expressed in 

terms of a unique rf- dimensional copula C as 

F{x) = C{fW(xi), . . . , FM(xd)}, X e W^. 



[Table 3 about here.] 



Table [3] reports rejection percentages under Hq when Fq is a bivariate c.d.f. with Clay- 
ton or Gumbel-Hougaard copula and standard exponential or standard normal margins. 
The parameter of the copula is chosen so that the theoretical value of Kendall's tau is 
equal to r e {0,0.25,0.5,0.75}. Note that r = corresponds to independence while 
r = 0.5 entails moderate dependence between the two components. Note also that the 
parameter m used in to discretize uniformly S2 was set to 8 (see Appendix [B] for more 
details). The settings m = 6 and m = 10 were also considered but this did not seem to 
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affect tlie results mucli. As tlie tests wliose approximate p-value is computed using tlie 
processes Dn ■* appeared systematically too liberal for the sample sizes under considera- 
tion, we only report the results of the tests based on the Dn^ in Table [31 As can be seen, 
the latter tests seem to hold their level reasonably well except perhaps the tests based 
on V^,v and which might be slightly too liberal for r = 0.75 and the sample sizes 
under consideration. Similar results were obtained in dimension three. 



[Table 4 about here.] 



Estimated rejection rates when both distributions in Hi have the same copula (Clay- 
ton or Gumbel-Hougaard) but differ in one margin are given in Table HI To be precise, 
the bivariate distributions Fi and F2 only differ in the first margin which is standard 
exponential for Fi and exponential with rate 0.5 for F2, while the second margin of both 
Fi and F2 is standard exponential. As one can see, the tests based on half-spaces are 
substantially more powerful than the tests based on lower-left orthants. The test based 
on f7„ V is the most powerful, overall, when r = 0, while, when r = 0.5, it is the test 
based on 14,+ that has, overall, the highest rejection rates. 



[Table 5 about here.] 



Table [5] reports rejection percentages when the change in distribution is only due to a 
change in the dependence structure: both Fi and F2 have standard exponential margins 
but the copula of Fi is the Clayton (resp. Gumbel-Hougaard) copula with a Kendall's tau 
of 0.1, while that of F2 is the Clayton (resp. Gumbel-Hougaard) copula with a Kendall's 
tau of T e {0.3,0.7}. As one can see from the overall low rejection percentages, this 
problem appears to be more difficult than the previous one. The tests have hardly any 
power for r = 0.3. For r = 0.7, the tests based on \4,v and 14,+ are the most powerful, 
overall, although one should be cautious as they might be slightly too liberal in the case 
of strongly dependent data according to Table O 

The setting used to obtain Table E] was finally extended to dimension three (results 
not reported). The conclusions are very similar to those obtained in dimension two with 
the difference that all the rejection rates are higher. Hence, as could have been expected, 
detecting a change in the dependence becomes easier as the dimension increases. 



5 Practical recommendations and illustration 



From the results of the Monte Carlo experiments partially reported in the previous sec- 
tion, the tests based on Sn,v and T„^+ seem good choices in the univariate case, while the 
test based on Vn,+ seems to be a good one in the multivariate case. 

As an illustration, we applied the tests based on the processes ©n^ t o the trivariate 



hydro logical data collected at the Ceppo Morelli dam, Italy, studied in ISalvadori et al. 



(120111 ). The data consist of annual maxima for 49 years (in the range 1937-1994) of 
three variables: L (dam reservoir water level in m), Q (peak flow in m^.s^^) and V (peak 
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vo lume in 10^ m^) . For a detailed description of the data, we refer the reader to Section 2 
of ISalvadori et all fcoilh . 



[Table 6 about here.] 

We first tested for a change in the distribution of (L,Q,V). The approximate p- values 
of the eight tests based on the processes are given in the first line of Table O Since 
there are both physical and stat i stical reasons to believe that L is independent of (Q,V) 
as explained in ISalvadori et al.l (120 111 ), as a second step, we tested for a change in the 
distribution of L and in the distribution of (Q,V) separately. The obtained approximate 
values are reported in the second and third lines of Table El As can be seen from the 
results of the test based on there is strong evidence of a change in the distributions 
of (L,Q,V) and (Q,V). From the second line of Table El we see that, on the contrary, 
there is no evidence of a change in the distribution of L. The latter finding is completely 
consistent with the fact that the variability of L is mainly due to the ma nagement policy 
of the reservoir which is constant over time. Indeed, as explained in ISalvadori et al. 



(120 111 ), the target of the dam manager is to keep a high water level in order to maximize 
electricity production. 

As classically done in the literature, under the hypothesis of a single break in a distri- 
bution, the change-point can be estimated by one of argmaxi^fcg;„_i Sn,k, argmaxi^fc^„_i T„ fc, 
argmaxi<jfc^„_i or eiigmaxi^k^n-i Vn,k depending on which test one wants to con- 
sider. For instance, the last estimator gives 31 for both {L,Q,V) and {Q,V), which 
corresponds to a change after the year 1976. 

Finally , let us mention that t he ap proach based on mul tivariate empirical c.d.f . s con- 
sidered in ICsorgo and HorvathI (Il997l . Section 2.6) and in iGombav and HorvathI (119991 ) 
has been extended by Inouel ( 2001 ) to serially dependent observations, although the lat- 
ter work is not aware of the former ones. A future research direction wou ld be to stud y 
generalizations as such proposed in this work in the setting considered by llnoud (120011 ) . 
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A Proofs 

A.l Proof of Theorem HI 

To prove Theorem [H we first show a lemma. 
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Lemma 1. Let be a P-Donsker class of functions with measurable envelope such 
that PF^ < 00. Then, 

sup - Z„(s,/)||j- ^ 0. 

se[0,l] 

Proof. For any s e [0, 1] and / e J^, we have 



Znis, f) ~~ Z„(s, /) = I ^ g 6 j (Pm/ ~ Pf)- 

Now, let Yfc = llPfc/ — -P/IIj^) ^ ^ {!)••• j''^}- Note that need not be measurable, 
but, by the assumption that sup^gjr < Fe{x) for all x e M'^, we have that 

W^'kf - Pf\W ^^kFe + PFe<^^^, 

where Pi-Fp + PFp is measurable. Thus, for each k, there exists a smallest random 



varia ble < oo such that Yk{(jj) ^ for every u (see Ivan der Vaart and Wellner 



2nnnl . Lemma 1.2.1). 

To prove the claim, it suffices to show that 



1 ^ 

An = max Y* X — V^i ^ 0, 



where the variables are all measurable. Now, for any n ^ 3, let = n ^/^(loglogri)^/^. 
Then, 

. '^k Ouk ^k Yf^ — 1 /9 1 1 "^A: 

An = max — X — — x — < max — x n ' log log n x max — , 

where 4 = ^"^I]f=i?*- 

Since J-" is P-Donsker with measurable envelope function F^ sat isfying PF^ < oo, we 



know from the law of the iterated logarithm for empirical processes (iDudley and Philipp 



19831 : iKosorokl . l2008l . page 31) that limsupy„*/a„ < oo almost surely, which implies that 
maxa^fe^fi Y'^*/afc ^ supj^^^Yj^ /ak < oo almost surely. Similarly, from the law of the 
iterated logarithm for the mean of an i.i.d. sequence with expectation and variance 1, 
we have max3^fc<c„ (^fc/a^ < ^'^Pk^s^k/o-k < °o almost surely. The desired result finally 
follows from the fact that n"^/^ log log n ^0. ■ 

Proof of Theorem [H. The proof of the weak convergence of the finite-dimensional 
marginal distributions of (Z„,Z„) to those of (Zp,Zp) is a more complicated version 
of the corresponding result for convergence of the rescaled random walk increments to 
Brownian motion. It is omitted here for the sake of brevity. 

To obtain that (Z„,Z„) (Zp,Zp) in {£°°([0,1] x J-")} ^, it remains to show that 



Z„,Z„) is asymptotically tight (see e.g. Ivan der Vaartl . Il998l . Section 18.3), which holds 



if both Z„ and Z„ converge weakly to tight random elements. 
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From Theorem 2.12.1 of Ivan der Vaart and Wellnerl (120001 ) . we have that Z„ 
X[0,1] X Now, let 



Zp in 



G' 



1 



be a multiplier version of G„. The class J-" being P- Donsker, from the functional un- 
conditional mu ltiplier central limit theorem (see e.g. Ivan der Vaart and Wellnerl . 12000 
Kosorokl . 2008 ). we have that G'p in i'^{J^), where G'p is an independent copy of 



the P-Brownian bridge Gp, which implies that G'^ is asymptotically tight. To show that 
Z„ Zp, where Zp is an independent copy of Zp, o ne can use the asymptotic tightnes s 
of and proceed as in the proof of Theorem 2.12.1 of lvan der Vaart and Wellner l2000h . 
Note that the pr oof can be further simplified if the process Z„ is symmetric in the sense 
of Chapter A.l of lvan der Vaart and Wellnerl (120001 ). which happens if ^i, . . . , are sym- 
metrically distributed aroun d zero. In that case, Ottav i ani's inequality can be replaced 
by Levy's inequality (see e.g. Ivan der Vaart and Wellnerl . l2000l Proposition A. 1.2), which 
shortens the argument. 

Hence, we have that (Z„,Z„) ^ (Zp,Z'p) in {£°°([0,1] x J^)}^. The fact that Zp is 
independent of Zp comes via that the finite-dimensional distributions, which are multi- 
variate normal, and the fact that Z„ and Z„ are uncorrelated. 



Next, notice that 



sup ||Z„(s,/) - Z„,(s,/)||j- 

se[0,l] 



IKf-Pfh^ sup 

se[0,l] 



[nsj 



* 1=1 



converges in outer probability to zero because sup^gj-Q j^j \n ' 21=1 s«l converges weakly 

to the supremum of the absolute value of Brownian motion and ||Pn/ — P/Wt ^ 0- 

From the continuous mapping theorem, we have that (Z„, Z„, Z„, Z„) -w (Zp, Zp, Zp, Zp) 
in {^°°([0,1] x J-")}^. The desired result finally follows from the previous remark and 
Lemma [H ■ 



A. 2 Proof of Theorem H 

Proof. For any s g [0, 1] and f e J^, let 

5^5, /) = {1 ~ A„(s)}Z„(s, /) - A„(s){Z„(l, /) - Z„(s, /)} = Z„(s, /) - Xn{s)Ml, /). 

Then, from Theorem [1], ([3]), the definitions of D„ and ©„, and the continuous mapping 
theorem, we obtain that 



D„,DJ ^ (©p„D'p^,Dy 
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in {£°°([0, 1] X J-")} . To obtain the desired result, it remains to show that Dn — '^n 
in £°°([0, 1] X J^). From the definitions of D„ and D„, we have that 



sup 

sg[0,1] 



K{s, f) -Bnis, f) < sup {1 - A„(s)} sup Z„(s, /) - Z„(s, /) 

^ se[0,l] se[0,l] 



+ sup A„(s) sup <{Z„(1,/) - Z„(s, /) - Z;(s,/) 

se[0,l] se[0,l] 



where Z* is defined in ([5]). The second supremum over s on the right of the previous 
inequahty converges outer almost surely to zero according to Lemma [H Furthermore, it 
can be verified that the last supremum over s (written for instance as a maximum over 
1 < \ns\ < n) is nothing else than the version of the second supremum computed from the 
"reversed" sequence (^^i^-^n), (Cn-ii -'^n-i), ■ ■ ■ ■, {ti^^i)- As the latter sequence and the 
original sequence have the same distribution. Lemma [1] implies that the last supremum 
converges in outer probability to zero. ■ 



A.3 Proof of Theorem [3] 

To prove Theorem [3l we first prove a lemma. 

Lemma 2. Let ^i, ■ ■ ■ be i.i.d. random variables with mean 0, variance 1 and satis- 
fying {P(|,^i| > x)}^^'^dx < 00, let Xi, . . . , Xn be i.i.d. random variables with law P 
independent of ^i, . . . ,C,n, o,nd let be a P-Donsker class with measurable envelope func- 
tion Ff. such that PFg < oo. Then, for every < ti < ^2 ^ 1 O'^d every e > 0, there exists 
M > such that 

[ns\ 

--VUiXd >M\<e. 
In '—^ 



supP* I sup 

nsN \sG[ti,t2] 



1=1 



Proof. Using the triangle inequality and the envelope Fg, for any s e [0, 1] and f s 
have 



[nsj 

1=1 



\ns\ 

' '^^ i=i 

Zn(-S,/) 



+ 



Yns\ 

In ■H^ 

i=i 



+ PFe X 



[ns\ 



'n ^ 

4=1 



Let < ti < t2 ^ 1- Since Z„, -w Zp in £°°([0, 1] x J^), we have that sup^gj-^^^^^] ll^n.(s, /)||jr 
is asymptotically tight. Let e > 0. Then, there exists M' > such that 







^ > M' j < e/4 


1 sup 




V se[ii,i2] 







It follows that, for n ^ 



1 sup 






V se[ti,t2] 
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Since sup^gj^^ j^] ll^n(s, f)\\jr is bounded by an almost surely finite random variable for any 
n < n^, there exists M" ^ M' such that 



1 sup 






VsG[ti,t2] 







Similarly, s n ^'^^Sllli' converges weakly to Brownian motion in £°°([0, 1]), and 
therefore sup^g^^^ j^] 1^^^"^^ Ij1=i' 61 is uniformly tight. Proceeding as above, there exists 
M'" > such that 



supP I sup 

wgN \se[ti,t2] 



^ 1=1 



> M'" <e/2. 



We get the desired result with M = M" + PF^ x M'". 

Proof of Theorem [3] (i). First, notice that, for any s g [0, 1] and any / e J^, 

Ktis,f) = 

Let K„,i(s) = ||A„(s) {1 - Us)} (P^„,j/ - ~ Kt{s, f)\\^. Clearly, 



sPJ s{tP^f + (1 - t)P^f} if s < t, 

tPl/ + (S - t)P2/ - s{tP,f + (1 - t)P2/} if S > t. 



sup Kn,t{s) =max< sup K„,t(s), sup K„,t(s) 

se[0,l] («6[0,t] se(t,l] 

Furthermore, for any s e [0, 1] and f ^ J^, 
Hence, 

sup K„,t(s) < sup \\Xn{s)F^nsif ~ + SUp ||A„(s)P„/ - s{tPif + (1 - t)P2/}||^. 

se[0,t] se[0,t] se[0,l] 

(7) 

The first term on the right of the previous inequality is smaller than 

sup ||A„(s) (P[nsj/ - Pif)\\jr + sup ||{A„(s) - s}Pi/||_^ 

se[0,t] se[0,t] 

= sup \\Zn{s, f)\\j, + sup |A„(s) ~ S| X 

V se[0,t] se[0,t] 

and therefore converges in outer probability to zero because sup^gp ^] ||Z„(s, /)||jr con- 
verges in distribution. The second supremum on the right of ([7j) is smaller than 



sup |A„(s) - S\ X ||P„/||^ + ||P„/ - {tPJ + (1 - t)P2f}\\ 

se[0,l] 



T 



and therefore converges in outer probability to zero because, T being Pi and P2-Donsker, 
P„/ converges in outer probability to tPi/ + (1 — t)P2/ uniformly in / e J-". 
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Similarly, 



sup K„,t(s) < sup ||A„(s)P^„,j/ - tPJ {s t)P2f\\^ 

se{t,l] se(t,l] 

+ sup \\Xnis)¥^f^s{tPJ + {l-t)P2f}y. 

se[0,l] 

We already know that the second supremum on the right converges to zero in outer prob- 
ability. Using the fact that, for s s {t, 1], A„(s)P[„,j = A„(t)P[„tj + {A„(s) - A„(t)}P[^^^']!^^^^j, 

where lP[„'sj!^|^„tj = {[ns\ — [^t-^J)^^ I]1=|Lj+i ^Xi, the first supremum is smaller than 



A„(t) {Pl„,j/ - PJ}\\^ + sup {A„(s) - A„(t)} {Fli:f\^^^f P2f] 

se{t,l] ^ ' 



T 



+ |A„(t) - i\ {\Pi]\j + + sup |A„(s) - s| X llPa/ll^ 

se(t,l] 

and converges to zero in outer probability because v^||A„(t){P[„tj/— Pi/}||jr = 

and V^sup^g(i 1] ||{A„(s) - A„(t)}{Pj';'j!J^^^j/ - P2/}i|j- converge in distribution. ■ 

Proof of Theorem [3] (ii) . Let us first show that sup^gp ij l|Z„(s, /)l|_7r is bounded in 
outer probability. Since 



sup 

sg[0,1] 



max < sup 

^ I SG[0,t] 





, sup 






^ SG(t,l] 





it is sufficient to verify the claim for each of the suprema in the maximum. From the 
definition of Z„, under Hi, the first supremum converges weakly from the continuous map- 
ping theorem and is therefore bounded in outer probability. For the second supremum, 
we write 



sup 

se(t,l\ 



Ms J) 



= max 

T [nt\<k^n 



n 



where 



and 



[nt\ 



Bn = max 



Cr, = max — — 



[nt\<k^n -v/n 



[nt\ 



j=[ntj+l 



Dn = max —= 

[nt\<k^n Jn 



( k \ . \nt\ 



/ j = [n,fj + l 
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The quantity An is clearly bounded in outer probability by Lemma [2] with ^1 = ^2 = t. 
For Bn, we can write 



Br, 



1 



max 



[nt 



J" 



< sup 

se[0,l-t] 



J" 



+ 



Ln(l-t)J 

^ 6+[?itj/(-^«+LntJ. 



J" 



where the two last terms on the right come from the fact that < [ntj — — t)J ^ 1. 
The first two terms on the right are bounded in outer probability by Lemma O and so is 
the third one because it converges almost surely to zero. It follows that B^ is bounded 
in outer probability. Now, C„ is bounded above by 



1 



max 



i=l 



\nt\ \ 

- y fix.) } max 



< max —j= 

[ntj<fc!£n -Jn 



[nt\ \ 



2& 



[nt\ 



X max 

[ntj<fc^n K 



where 



= sup 

se[0,l] 



1 



Similarly, we have that D„ < /„. To show that C„ and D„ are bounded in outer proba- 
bility, we will show that /„ is bounded in probability. 

Since Fe{Xj) converges to tPiFe + (1 — t)-P2-Fe < °o almost surely, we have 

that sup„gi^ 2j=i ^e{Xj) is an almost surely finite random variable. The fact that 
is bounded in probability then follows from the weak convergence of s t-^ n^^^'^ Xj1=i' 
to Brownian motion in £'^([0, 1]), which implies that sup^g^ 1^"^''^ I]l=i' 61 is uniformly 
tight. 

Thus, An, Bn, Cn, and D„ are all bounded in outer probability, which implies that 
suPsg(t^i] ||Zn('S, /)|| J- is bounded in outer probability, and therefore that sup^g^ ||Zn(s, 
is bounded in outer probability. 

The analogous result for the process Z* follows from the fact that sup^g^ W^ni^, 
(written for instance as a maximum over 1 < [ns\ < n) is nothing else than the version 
of sup^g[oi] ||Z„(s, computed from the sequence {^n,Xn),{^n-i,Xn-i), . . . ,{Ci,Xi) 
which has the same distribution as the original sequence. The desired result is finally an 
immediate consequence of (jl]). ■ 



Proof of Theorem |3] (iii). For any s e [0, t] and / e J^, we can write 

[ns\ [ns\ 
Zn{s, /) = ^ 2 UfiXd Pi/} - {Fnf Pi/} X ^ 2 

* 1=1 * 4=1 
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while, for any s e [t, 1] and / g J^, we have 



[nt] [nt\ 



[nsj [nsj 



j=[ntj + l ^ i=[ntj + l 



Now, let us show that 



( \ns\ \ns\ \ 

(s, /) - ^ 2 ?.{/(^.) - A/}, ^ 2 ^0 ^ ^) ^ ^i(^)) 

* 1=1 * 1=1 / 

in {£'*([0,t] X J-")}^, where i?i is Brownian motion and is independent of Zp^. Using the 
multivariate central limit theorem, it can be verified that we have weak convergence of 
the finite-dimensional distributions. Joint asymptotic tightness follows from the weak 
convergence of s > rr^l'^'^^^\ to Bi in £^^([0, 1]), and the weak convergence of Z„ to 
Zp^ in £°°([0, t] X J^). Since the two component processes on the left of Qj are uncorrelated, 
their weak limits are independent. Similarly, we have that 



( [ns\ [ns\ \ 

(sj)^^ Uf{x,)-p,f},{s,f)^— J] 
^ i=[n.ij + l ^ i=[nt\ + l / 



^ ((s, /) - Zp,(s - /), {s, f) ^ B,{s ~ t)) (10) 

in 1] X J-")}^, where B2 is Brownian motion independent of Zpj, and (Zp2,i?2) is 

independent of (Zpj,i?i). 

Combining the fact that P„,/ converges in outer probability to tPif + (1 — t)P2f 
uniformly in / e with ([9]), we obtain from the continuous mapping theorem that, in 
e^i[0,t] X Z„ converges weakly to (sj) ^ Wp,(s,/) = Zp,(s,/) - (1 - t){P2f 
Pif)Bi{s). Using (ini similarly, it can be verified that, in ^^([t, 1] x J^), 

[nsJ [nsJ 



n ■\/n , , 

i=[nij + l ^ i=[nt\ + l 

converges weakly to (s, /) ^ Wp,(s,/) = Zp,(s - t, /) - t{P2f Pif)B^{s t). By 
independence of first two terms in ([HD with the last two terms, we then obtain that, in 
1] X J^), Z„ converges weakly to (s, /) ^ Wp^{t, f) + Wp^is, /), which implies that 
Z„ converges weakly in £°°([0,1] x J^). The desired result finally follows from the fact 
that ©ri(s, /) = Z„(s, /) — A„(s)Z„(l, /) and the continuous mapping theorem. ■ 

A. 4 Proofs of Propositions [1] and [2] 

We state a lemma before giving the proofs of the propositions. 
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Let A be the space of bounded Borel measurable functions on M'^ and let B be the 
space of c.d.f.s of Borel probability measures on M.'^. The spaces A and B are subsets 
of £°°(M ) and the topologies on A and B are the ones induced by uniform convergence. 
The following result (the help of Johan Segers is gratefully acknowledged) will allow us 
to apply the continuous mapping theorem in the proof of Proposition [T] 

Lemma 3. Let (p : A x B ^ M. be defined by (j){a, b) = S^^^ o-db. The map (p is continuous 
at each {a,b) e A x B such that a is continuous on Mf^. 

Proof. Let (a„,6„) be a sequence in ^4 x S such that sup^jgi^d |a„(x) — a{x)\ and 
sup^jgigd \bn{x) — b{x)\ 0. It is sufficient to show that ^jgda^db^ ^jg^adfe. By the 



On dbn ~ a dbn + a dbn — adb 

jRd Jr<* jRd 



triangle inequality, 

On d6„ - a db 
For the first term on the right of this inequality, we have 

andbn- adbn ^ \ |a„ - a| dfe^ 

jRd jRd jRd 

< sup \an{x) — a{x)\ dh, 



sup \an[x) — a[x)\ 0. 



For the second term, since sup2,gigd \l)n{x) — b{x)\ 0, we have that bn{x) —> b{x) for every 
X G M"', which, by the Portmanteau lemma and the continuity of the function a implies 
that a dbn ^r^ ^ ^6. ■ 

Proof of Proposition [T]. A consequence of Theorem [2] and the fact that, under Hq, F„ 
converges almost surely to Fq (the c.d.f. of Pq) uniformly in a; e M"^ is that 



1 ) • • • ) ' i ■ ■ ■ i ' ^ n 



in {£°°([0, 1] X ]R'^)}(^^+^). Using the map (f) defined in Lemma [31 it is easy to see 
that Sn,v = ^^Pse[o,i]4>[{^n{s,-)y,Fn], that Sii^l = sup^g[o,i] ^[{^"^(s, Ol^'-^nJ and that 
Sn^v = supj,g|-Q 0[{Dn ''(s, ■)}^, F„] , j G {1, . . . , A^}. Furthermore, the limiting process 
is continuous almost surely. The result then follows from Lemma [3] and the continuous 
mapping theorem. ■ 

Proof of Proposition [2l From Theorem |3l we have that 



sup sup 

se[0,l] xeW^ 



i{s,x) 



n 



Kt{s,x) 



where Kt{s, x) = {Fi{x) - F2{x)}{s a - (s v t)}, i.e., D^/V^ ^ Kt in ^"([0, 1] x M ). 
Also, under Hi, Fn converges almost surely to Ft = tFi + (1 — t)F2 uniformly in a; 6 M"'. 
Then, 

{On(s,x)}2 



sup 

sg[0,1] 



Ir 



n 



-dFJx) 



\ {Kt{s,x)fdFt{x) 
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^ An + Bn, 



11) 



where 



and 



We have 



An = sup 
se[0,l] 



Bn = sup 

sg[0,1] 



n 



dFJx) 



dFJx] 



\ 



An < sup sup 

sg[0,1] x-eIRd 



{Kt{s,x)YdFn{x)- [ {Kt{s,x)YdFt{x) 

{Dn(s,x)}2 



n 



{Kt{s,x)y 







by the continuous mapping theorem, and, with the notation Qs = {Kt{s, ■)}^, s s [0, 1], 

Bn = sup \Xnit)]P[nt\gs + {1 - Xnit)}K-[nt\9s ~ tPiQs - (1 - t)P2^7s}| 
sg[0,1] 

\n{t) sup \{f>ynt\ - Pl)9s\ + {1 - A„(t)} SUp \{K-[nt\ ' ^2)^s| + 2|A„(t) - t| 



se[0,l] 



sg[0,1] 



because sup^g^ ^ | Pig's] ^ 1 and sup^gp^ij \P2gs\ ^ 1- Now, 

' ' nt\ - Pi)gs\ ^ |(Pm - Pi){Fr F2f\ ^ 



sup 

as[0,l] 



and 



sup |(P:_^„,j ~ P2)gs\ < |(P:-^„*j - P2)(Pi - F,f\ ^ 

se[0,l] 



by the law of large numbers, which implies that S„ -> 0. It follows from (1111) that 

{^n{s,x)Y 



-dFn{,x) 



{Kt{s,x)YdFt{x) 



in £°°([0, 1]), from which, using the continuous mapping theorem, we obtain that 



5'n,v/?^ sup {Kt{s,x)] dFt{x). 

ss[0,l] JK.'i 

Since Fi and F2 are distinct and Kt{s, x) = {Fi{x) — F2{x)]{s a t){l — (s v t)}, we have 
that sup^gp^i] S^^aiKtis, x)YdFt{x) > 0, which implies that Sn,v ^ +00. 

Now, let j G {1, . . . , A^}. The claim for Sn^v follows from the inequality 



5^^^= sup f 0^Xs,x)ydFn{x)^ sup sup{Di^)(s,y)}' f 

sg[0,1] JM.d' se[0,l] yeRd. J^, 



dFnix) < 



and the fact that Tn,v is bounded in outer probability from Theorem [31 It remains 



finally to prove the claim for Sn^v A consequence of Theorem [3] is that (Dn',Fn 



converges weakly in {£°°([0,1] x M"^)}^, which, combined with the fact that Sn}~/ = 
suPsg[o,i] 0[{Dn ^s, •)}^, F„] and Lemma O implies that Sn}v converges in distribution, 
and hence, that it is bounded in outer probability. ■ 
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B Implementation of the tests based on Un,k cind T4 



k 



Let m > 2 be an integer and let be elements of uniformly spaced over 

For k E {1, . . . ,n — 1}, the following numerical approximations are then considered: 



Un,k^ -y,-y,hn{-,ai,aixA 

1=1 9=1 ^ ^ ^ ^ 

1 = 1 q=l y 4=1 i=fc + l 



and 



Vnk ^ max max 

ligi^m l^q^n 



A;(n — A;) 

-7- — max max 



j=i 



n — k 



Next, recall that, for any j e {1, . . . , A^}, (resp. ©n^) is the version of D„ (resp. D„) 
computed from \ . . . , Proceeding as above, for any j e {1, . . . , A^}, the multiplier 
versions of Un,k and Vn^k based on the process Dn ■* are computed respectively as 



U 



n,k 



-, m n 

J_yy 



m n I k 

[n — 



1=1 q=l 



j=l 



n 



2 



and 



^ n '^''^ max max 

liSisSm IsSgsCn 



n 



j=fc+l 



i=l 



i=fc+l 



Similarly, for any j s {1, . . . , A^}, the multiplier versions of Un,k and Vn,k based on the 
process Dn ■* are computed respectively as 



U 



n,k 



_J_yy 



1=1 g=l 



[n — 



j=i 



n 

-k 2 i^\l{alX,^aJX,)-Fa,A^JX,)] 



i=k+\ 
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and 



k 




3/2 



max max 



n 



j=fc+l 



where „ is the empirical c.d.f. computed from the projected sample a Xi, . . . , a Xn- 
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Table 1: Percentage of rejection of Hq in the univariate case computed from 1000 random 
samples of size n generated under Hq defined in ([1]) where Fq is the c.d.f. of the standard 



normal distribution. 



n 






c* 

n, V 








T 


T 


n, V 


Tn,+ 


Tn,+ 




50 


7.2 


5.7 


5.7 


7.7 


5.1 


5.9 


6.7 


5.8 


5.7 


8.4 


5.2 


4.3 


100 


6.5 


5.5 


6.2 


6.0 


4.9 


6.1 


7.1 


6.6 


6.5 


8.1 


6.2 


6.2 


200 


5.9 


4.8 


4.3 


5.9 


4.6 


5.5 


5.4 


4.5 


4.5 


7.5 


5.4 


4.3 



25 



Table 2: Percentage of rejection of Hq in the univariate case computed from 1000 samples of size n generated under Hi defined 
in ([2]), where k* = [nt\, and Fi and F2 are the c.d.f.s of the distributions given in the first two columns. 





F2 




n 




t 






*^n, V 


'S'n,+ 


Sn,-\- 






-"^ n, V 




V 


Tn,+ 


Fn,+ 


Fn,+ 


N 


;o,i) 


N 


;o.5,i) 


50 





10 


9.1 


7.1 


6.8 


10.1 


7.1 


8.2 


9.0 


7.3 


6 


9 


14.2 


8.3 


8.6 


N 


;o,i) 


N 


;o.5,i) 


50 





25 


18.4 


15.9 


13.4 


19.1 


14.4 


14.4 


17.0 


13.6 


8 


9 


23.0 


13.6 


12.0 


N 


;o,i) 


N 


;o.5,i) 


50 





50 


34.0 


30.5 


31.6 


34.1 


29.0 


32.1 


32.6 


29.5 


25 


6 


35.4 


26.9 


24.9 


N 


;o,i) 


N 


;o.5,i) 


100 





10 


9.7 


8.2 


6.5 


10.1 


8.5 


7.3 


9.7 


8.7 


6 


7 


11.6 


9.4 


7.3 


N 


;o,i) 


N 


;o.5,i) 


100 





25 


36.9 


34.0 


33.8 


36.2 


33.3 


31.7 


30.4 


28.6 


25 


1 


37.2 


31.4 


29.3 


N 


;o,i) 


N 


;o.5,i) 


100 





50 


58.7 


55.9 


54.1 


56.7 


53.9 


56.0 


52.1 


48.6 


42 


1 


53.6 


49.3 


45.5 


N 


;o,i) 


N 


;o.5,i) 


200 





10 


15.9 


15.1 


14.0 


18.6 


17.3 


16.8 


15.3 


14.6 


11 


1 


22.1 


19.5 


18.0 


N 


;o,i) 


N 


;o.5,i) 


200 





25 


65.2 


64.1 


63.9 


65.8 


64.3 


62.5 


56.7 


55.5 


54 


8 


61.6 


57.2 


52.2 


N 


;o,i) 


N 


;o.5,i) 


200 





50 


87.1 


86.4 


86.8 


86.1 


85.3 


85.0 


81.2 


80.6 


79 


1 


81.2 


79.5 


76.2 


N 


;o,i) 


N 


;o,2) 


50 





10 


6.8 


5.6 


5.9 


8.1 


5.4 


4.9 


8.1 


6.0 


5 


1 


11.9 


7.0 


6.3 


N 


;o,i) 


N 


;o,2) 


50 





25 


8.9 


6.4 


7.1 


11.7 


9.4 


6.9 


11.3 


10.0 


8 


1 


17.0 


10.9 


8.4 


N 


;o,i) 


N 


;o,2) 


50 





50 


12.9 


10.3 


9.9 


18.8 


13.5 


14.1 


19.6 


16.6 


12 


6 


27.8 


18.4 


17.4 


N 


;o,i) 


N 


;o,2) 


100 





10 


7.1 


6.4 


6.6 


8.1 


6.8 


6.3 


8.6 


7.7 


6 





10.7 


8.0 


6.0 


N 


;o,i) 


N 


;o,2) 


100 





25 


8.8 


7.5 


8.1 


17.5 


14.3 


15.4 


14.2 


13.2 


12 


1 


25.3 


19.2 


18.8 


N 


;o,i) 


N 


;o,2) 


100 





50 


23.1 


20.8 


24.9 


38.4 


34.2 


39.1 


34.5 


31.9 


27.5 


46.6 


38.9 


37.7 


N 


;o,i) 


N 


;o,2) 


200 





10 


6.4 


5.9 


5.9 


9.1 


8.2 


6.9 


9.0 


8.7 


6 


7 


14.5 


12.2 


8.6 


N 


;o,i) 


N 


;o,2) 


200 





25 


21.0 


19.6 


21.2 


50.2 


47.8 


41.3 


33.6 


32.6 


32 


5 


54.7 


49.3 


43.8 


N 


;o,i) 


N 


;o,2) 


200 





50 


64.8 


64.0 


65.1 


82.8 


81.7 


82.3 


72.1 


71.0 


65 


7 


83.7 


82.5 


79.1 


E( 


1) 


E( 


0.5) 


50 





10 


8.1 


6.8 


6.0 


9.3 


7.2 


7.3 


8.9 


7.4 


5 


8 


14.5 


8.7 


7.8 


E( 


1) 


E( 


0.5) 


50 





25 


28.0 


25.3 


22.1 


30.7 


25.5 


22.6 


25.5 


22.4 


19 


6 


33.8 


24.2 


22.3 


E( 


1) 


E( 


0.5) 


50 





50 


50.5 


46.2 


46.3 


49.6 


43.8 


48.8 


45.0 


40.8 


40 


7 


48.6 


40.7 


42.9 


E( 


1) 


E( 


0.5) 


100 





10 


12.7 


11.0 


7.9 


14.3 


11.7 


11.4 


12.8 


11.6 


11 


1 


16.7 


13.1 


10.9 


E( 


1) 


E( 


0.5) 


100 





25 


51.9 


50.0 


50.8 


52.2 


48.4 


47.4 


45.2 


42.9 


41 


3 


49.9 


44.9 


45.3 


E( 


1) 


E( 


0.5) 


100 





50 


77.4 


76.1 


74.6 


77.3 


74.8 


72.9 


72.4 


70.4 


70 


7 


73.8 


69.4 


67.9 


E( 


1) 


E( 


0.5) 


200 





10 


19.9 


19.3 


20.2 


24.9 


23.2 


23.8 


16.8 


16.3 


13 


9 


25.4 


23.3 


19.4 


E( 


1) 


E( 


0.5) 


200 





25 


87.0 


86.2 


85.5 


84.4 


83.1 


81.4 


79.0 


78.0 


76 


2 


80.8 


77.3 


77.9 


E( 


1 


E( 


0.5) 


200 





50 


96.6 


96.3 


95.7 


95.4 


94.9 


95.7 


94.9 


95.0 


94 


7 


93.2 


92.7 


91.7 



Table 3: Percentage of rejection of Hq computed from 1000 random samples of size n generated under Ho defined in ([T]), where Fq 
is a bivariate c.d.f. whose univariate margins F^^^ and Fq^^ are either both standard normal (N) or both standard exponential (E), 
and whose copula is either the Clayton (CI) or the Gumbel-Hougaard (GH) with a Kendall's tau of r. The parameter m used to 
uniformly discretize S2 was set to 8 (see Appendix [B] for more details). 

Ci GH 





n 




r 






T 




Un,v 


Un,+ 




Vn,+ 


Sn, V 




T 




Uji, V 


Un,+ 




Vn,+ 


N 


50 


0, 


.00 


3.8 


3.7 


4.7 


4.8 


5.5 


4.4 


5.9 


5.9 


4.2 


3.8 


5.4 


4.0 


4.7 


3.5 


5.8 


5.7 


N 


50 


0, 


.25 


3.2 


2.6 


4.8 


4.0 


5.5 


4.1 


6.3 


5.7 


4.9 


4.4 


4.5 


4.2 


5.2 


4.4 


6.2 


5.4 


N 


50 


0, 


.50 


5.3 


4.9 


6.3 


6.2 


4.9 


4.4 


5.1 


5.4 


4.7 


4.4 


5.8 


5.3 


4.9 


4.3 


5.3 


5.2 


N 


50 


0, 


.75 


5.0 


4.7 


5.0 


4.6 


4.9 


4.5 


6.4 


5.7 


4.4 


3.9 


4.5 


4.1 


4.2 


3.8 


4.6 


4.5 


N 


100 


0, 


.00 


4.9 


4.9 


5.5 


6.2 


5.6 


4.7 


6.5 


6.6 


4.0 


4.5 


5.2 


4.7 


3.5 


3.4 


4.0 


4.2 


N 


100 


0, 


.25 


5.4 


4.9 


5.6 


5.1 


4.7 


4.3 


5.6 


5.6 


6.0 


4.9 


5.9 


6.3 


5.3 


4.4 


6.5 


6.1 


N 


100 


0, 


.50 


5.1 


4.2 


5.1 


4.9 


4.1 


4.5 


6.0 


6.3 


4.7 


4.2 


4.3 


3.9 


3.8 


4.2 


5.7 


5.4 


N 


100 


0, 


.75 


4.9 


5.4 


5.5 


5.7 


5.1 


5.3 


6.5 


7.4 


5.1 


4.7 


4.5 


3.6 


5.5 


4.2 


5.3 


5.1 


N 


200 


0, 


.00 


5.5 


5.4 


5.5 


6.0 


5.1 


6.0 


5.7 


7.4 


4.8 


4.2 


5.6 


4.6 


4.0 


3.7 


5.3 


4.7 


N 


200 


0, 


.25 


5.4 


5.2 


5.9 


6.4 


4.9 


4.4 


6.0 


6.3 


5.2 


4.1 


5.3 


5.2 


4.9 


4.3 


5.7 


4.8 


N 


200 


0, 


.50 


5.8 


4.8 


4.9 


4.1 


5.2 


4.4 


4.9 


4.9 


5.1 


6.0 


5.5 


6.0 


6.0 


6.0 


5.8 


7.1 


N 


200 


0, 


.75 


6.7 


5.5 


7.3 


6.0 


6.7 


6.1 


7.1 


6.3 


4.9 


4.5 


5.4 


4.7 


5.8 


4.6 


5.6 


5.6 


E 


50 


0, 


.00 


4.3 


4.8 


6.0 


6.0 


4.0 


3.9 


6.1 


6.1 


3.8 


3.9 


3.8 


3.6 


3.6 


3.2 


5.8 


5.7 


E 


50 


0, 


.25 


3.7 


3.4 


5.3 


4.0 


4.7 


3.6 


5.8 


4.4 


6.6 


5.6 


6.9 


5.9 


5.2 


4.2 


6.7 


5.8 


E 


50 


0, 


.50 


3.7 


3.4 


4.8 


5.0 


3.3 


3.4 


6.5 


6.2 


4.8 


4.7 


5.9 


5.1 


4.7 


3.6 


6.9 


5.4 


E 


50 


0, 


.75 


5.4 


4.9 


6.0 


5.4 


5.2 


4.2 


5.9 


6.4 


5.9 


5.2 


6.0 


5.2 


6.0 


5.2 


7.6 


6.2 


E 


100 


0, 


.00 


4.9 


4.6 


5.3 


4.6 


4.6 


4.8 


4.7 


5.1 


5.3 


4.5 


5.2 


4.6 


5.5 


4.6 


5.5 


5.7 


E 


100 


0, 


.25 


5.1 


4.9 


4.9 


4.7 


5.9 


4.8 


5.7 


6.1 


3.5 


3.6 


4.4 


4.6 


4.1 


4.6 


5.6 


6.5 


E 


100 


0, 


.50 


5.4 


4.6 


5.3 


4.5 


5.3 


4.4 


6.3 


5.1 


5.2 


4.4 


4.0 


3.8 


4.9 


3.8 


5.5 


5.9 


E 


100 


0, 


.75 


5.5 


5.8 


5.5 


5.6 


5.1 


5.0 


5.3 


5.6 


6.1 


6.5 


7.2 


7.1 


6.6 


6.0 


6.9 


6.8 


E 


200 


0, 


.00 


4.9 


4.4 


5.0 


4.9 


5.2 


4.6 


5.1 


5.4 


5.1 


4.4 


4.6 


4.9 


4.5 


4.7 


5.5 


5.6 


E 


200 


0, 


.25 


5.6 


5.6 


5.4 


5.3 


5.6 


5.3 


6.6 


6.4 


4.5 


4.6 


5.5 


5.0 


4.6 


4.9 


4.7 


5.7 


E 


200 


0, 


.50 


5.5 


4.8 


6.4 


5.3 


5.1 


4.8 


5.8 


5.5 


6.7 


6.2 


6.3 


6.0 


6.1 


6.2 


5.6 


6.5 


E 


200 


0, 


.75 


4.9 


4.3 


4.9 


5.0 


4.7 


4.6 


6.2 


6.3 


6.8 


5.6 


5.8 


5.5 


6.7 


5.9 


7.4 


6.7 



Table 4: Percentage of rejection of Hq computed from 1000 samples of size n generated under Hi defined in ([2]), where k* = [nt\, 
Fi and F2 are bivariate c.d.f.s that only differ in their first margin which is standard exponential for Fi and exponential with rate 
0.5 for F2. The common copula C of Fi and F2 is either the Clayton (CI) or the Gumbel-Hougaard (GH) with a Kendall's tau of 
r. The second margin of both Fi and F2 is standard exponential. The parameter m used to uniformly discretize S2 was set to 8. 

Ci GB. 



n 




r 




t 


Sn, V 




T 










Vn,+ 


>~'n, V 




T 

-'■ n, V 










K,+ 


50 





.0 


0, 


.10 


5.3 


5.1 


5.7 


5.4 


6.1 


6.8 


7.4 


9.0 


5.6 


5.6 


6.5 


7.2 


6.8 


7.2 


7.9 


9.0 


50 





.0 


0, 


.25 


12.4 


12.4 


12.3 


13.8 


19.6 


19.1 


16.6 


18.4 


10.6 


11.0 


12.8 


13.4 


19.4 


19.5 


17.4 


19.1 


50 





.0 


0, 


.50 


21.3 


19.9 


23.8 


22.2 


41.2 


36.7 


35.0 


35.2 


24.6 


23.4 


25.2 


23.5 


43.4 


39.5 


35.2 


35.4 


50 





.5 


0, 


.10 


4.3 


3.9 


4.8 


4.8 


4.3 


4.8 


6.4 


8.4 


5.0 


5.2 


5.9 


5.5 


7.3 


7.7 


9.3 


10.6 


50 





.5 


0, 


.25 


8.3 


7.9 


11.6 


12.1 


19.5 


20.8 


26.3 


30.2 


8.7 


8.7 


11.8 


11.8 


21.6 


25.2 


30.0 


32.9 


50 





.5 


0, 


.50 


17.5 


15.3 


26.1 


22.9 


49.2 


50.1 


59.5 


56.8 


18.6 


18.2 


26.9 


23.9 


48.5 


48.7 


63.1 


58.1 


100 





.0 


0, 


.10 


7.9 


8.3 


7.9 


9.7 


9.8 


10.4 


9.9 


11.9 


7.3 


7.0 


8.7 


9.3 


10.1 


10.6 


9.1 


11.6 


100 





.0 


0, 


.25 


25.0 


23.9 


23.7 


26.1 


44.5 


42.7 


33.5 


35.9 


25.9 


25.9 


27.5 


30.6 


47.3 


45.3 


34.6 


40.0 


100 





.0 


0, 


.50 


44.1 


43.1 


49.9 


47.7 


72.1 


68.6 


63.8 


61.6 


45.8 


44.6 


53.5 


49.8 


74.2 


70.0 


63.8 


62.3 


100 





.5 


0, 


.10 


5.9 


5.5 


6.4 


7.3 


8.9 


9.7 


11.6 


16.4 


6.6 


6.6 


6.9 


8.4 


9.3 


11.1 


10.8 


18.4 


100 





.5 


0, 


.25 


14.7 


15.6 


22.4 


24.5 


50.8 


53.7 


57.5 


60.3 


17.9 


17.3 


28.1 


28.0 


59.8 


61.6 


70.1 


70.5 


100 





.5 


0, 


.50 


28.2 


30.2 


49.0 


48.7 


84.7 


85.8 


89.2 


87.6 


34.2 


34.3 


51.4 


48.6 


87.7 


88.8 


94.2 


92.6 


200 





.0 


0, 


.10 


11.2 


11.5 


10.5 


14.4 


15.7 


19.1 


11.6 


19.3 


11.2 


12.0 


10.3 


13.9 


17.4 


19.7 


12.8 


19.4 


200 





.0 


0, 


.25 


49.6 


48.4 


56.9 


57.1 


81.2 


77.5 


67.5 


68.6 


52.6 


51.5 


58.6 


58.4 


81.4 


79.0 


67.0 


71.0 


200 





.0 


0, 


.50 


80.1 


78.5 


87.4 


84.9 


97.0 


95.8 


92.6 


91.6 


77.3 


75.9 


85.9 


81.9 


97.2 


96.2 


92.1 


90.4 


200 





.5 


0, 


.10 


7.0 


7.8 


9.9 


14.6 


16.3 


22.9 


21.1 


30.2 


8.0 


8.4 


8.7 


12.3 


18.6 


24.3 


22.2 


35.2 


200 





.5 


0, 


.25 


25.8 


28.3 


53.6 


54.2 


87.7 


88.7 


90.1 


90.4 


30.4 


32.8 


55.6 


54.7 


95.0 


94.6 


97.7 


95.9 


200 





.5 


0, 


.50 


52.9 


57.0 


83.9 


82.3 


99.4 


99.3 


99.4 


99.5 


67.5 


66.9 


86.7 


84.7 


99.9 


99.8 


100.0 


100.0 



Table 5: Percentage of rejection of Hq computed from 1000 samples of size n generated under Hi defined in (|2]), where k* = [nt\, Fi 
and F2 are bivariate c.d.f.s with standard exponential margins and copula from the same family (Clayton (CI) or Gumbel-Hougaard 
(GH)) but with different parameter values. The copula of Fi has a Kendall's tau of 0.1, while that of F2 has a Kendall's tau of r. 
The parameter m used to uniformly discretize S2 was set to 8. 

Ci GB. 



n 




r 




t 


Sn, V 




T 










K,+ 






T 

n, V 


Tn,+ 








K,+ 


50 


0, 


.3 


0, 


.10 


5.5 


5.9 


5.2 


5.1 


4.3 


3.8 


5.2 


6.0 


4.7 


5.2 


5.1 


5.1 


5.8 


4.9 


6.8 


5.8 


50 


0, 


.3 


0, 


.25 


5.9 


5.1 


6.8 


5.3 


5.4 


4.7 


6.2 


6.5 


5.6 


5.1 


6.0 


5.6 


5.9 


4.2 


7.2 


6.2 


50 


0, 


.3 


0, 


.50 


4.4 


4.5 


6.8 


4.6 


3.9 


2.9 


6.2 


4.9 


6.1 


6.1 


6.3 


6.3 


4.7 


4.7 


6.8 


6.7 


50 


0, 


.7 


0, 


.10 


5.7 


5.7 


6.3 


6.2 


5.1 


4.9 


6.4 


6.4 


5.4 


5.6 


5.8 


5.6 


3.8 


4.2 


6.3 


6.6 


50 


0, 


.7 


0, 


.25 


7.5 


8.0 


6.8 


8.1 


4.4 


3.8 


6.9 


9.3 


7.3 


9.5 


9.1 


8.9 


5.8 


5.5 


8.0 


9.8 


50 


0, 


.7 


0, 


.50 


13.4 


13.7 


13.9 


14.3 


7.1 


8.3 


13.9 


14.8 


12.8 


13.1 


12.6 


11.8 


6.1 


6.9 


13.3 


14.0 


100 


0, 


.3 


0, 


.10 


4.6 


4.8 


5.9 


5.8 


4.4 


4.5 


6.5 


6.4 


5.0 


5.2 


5.9 


6.4 


6.4 


5.7 


8.0 


7.6 


100 


0, 


.3 


0, 


.25 


4.8 


5.1 


5.4 


5.3 


4.7 


4.4 


5.2 


5.8 


5.3 


5.6 


6.4 


7.2 


4.8 


4.3 


5.4 


6.4 


100 


0, 


.3 


0, 


.50 


7.6 


7.4 


7.2 


6.5 


5.3 


4.4 


6.2 


6.1 


7.0 


6.8 


7.0 


7.0 


4.4 


4.5 


5.8 


6.9 


100 


0, 


.7 


0, 


.10 


7.1 


6.5 


7.2 


6.8 


5.6 


4.8 


5.8 


8.7 


5.1 


5.8 


5.9 


5.8 


3.8 


5.0 


5.4 


7.0 


100 


0, 


.7 


0, 


.25 


12.9 


13.4 


13.0 


13.4 


6.9 


8.5 


12.0 


18.6 


13.9 


15.7 


12.7 


13.6 


5.8 


7.7 


13.5 


21.2 


100 


0, 


.7 


0, 


.50 


25.4 


26.7 


23.8 


25.0 


11.4 


15.0 


32.5 


36.0 


27.1 


27.1 


25.0 


25.7 


12.7 


16.6 


34.2 


39.4 


200 


0, 


.3 


0, 


.10 


6.0 


5.9 


5.0 


5.2 


5.6 


4.5 


5.4 


4.6 


4.5 


4.4 


4.9 


5.3 


4.9 


4.9 


6.0 


6.2 


200 


0, 


.3 


0, 


.25 


6.2 


7.1 


5.9 


6.9 


4.7 


5.6 


5.7 


7.1 


6.0 


6.2 


5.8 


5.5 


4.7 


5.3 


6.2 


6.8 


200 


0, 


.3 


0, 


.50 


7.9 


8.4 


7.4 


7.4 


5.9 


5.5 


6.9 


7.2 


10.8 


10.6 


9.3 


10.4 


7.2 


7.2 


8.6 


9.3 


200 


0, 


.7 


0, 


.10 


8.9 


11.2 


9.1 


11.5 


6.5 


7.1 


8.0 


10.8 


8.9 


9.6 


8.7 


11.1 


5.1 


6.7 


8.1 


11.8 


200 


0, 


.7 


0, 


.25 


26.0 


28.1 


25.3 


25.3 


11.0 


15.0 


25.1 


29.6 


28.6 


30.0 


25.9 


28.3 


14.3 


20.4 


30.1 


43.0 


200 


0, 


.7 


0, 


.50 


43.5 


44.8 


43.7 


44.4 


22.6 


30.4 


56.8 


61.4 


48.0 


45.9 


45.2 


44.5 


35.2 


51.7 


79.8 


85.2 



Table 6: Approximate p-values of the tests based on the processes Dn for the trivariate 
hydrological data considered in Section [51 The trivariate (resp. bivariate) tests based 
on half-spaces were run with m = 32 (resp. m = 8). The approximate values were 
computed from = 10, 000 multiplier realizations. 





Variables Sn,v 


\^ 


T 












(L,Q,V) 0.114 


0.120 


0.015 


0.028 


0.010 


0.010 


0.004 


0.006 


L 0.479 


0.314 


0.510 


0.236 










(Q,V) 0.024 


0.028 


0.012 


0.012 


0.015 


0.014 


0.004 


0.007 
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